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Abstract 



^^ , We prove tight bounds for the oo-norm of the inverse of a symmetric, diagonally dominant positive 

matrix J; in particular, we show that ||J _1 |[oo is uniquely maximized among all such J. We also prove 
a new lower-bound form of Hadamard's inequality for the determinant of diagonally dominant positive 

i-Q matrices and an improved upper bound for diagonally balanced positive matrices. Applications of our 

results include numerical stability for linear systems, bounds on inverses of differentiable functions, and 
consistency of the maximum likelihood equations for random graph distributions. 

1 Introduction 
> 

CN ' An n x n real matrix J is diagonally dominant if 

3§; A ? (J):HJ„;|-^|J ?;j |>0, for* 

m 

Irreducible, diagonally dominant matrices are always invertible, and such matrices arise often in theory and 
applications. As a recent example, the work of Spiclman and Teng [T31 Q3] shows how to solve diagonally 
dominant linear systems in nearly-linear time in the input size, a fundamental advance in algorithmic com- 
plexity theory and numerical computation. In this paper, we study the mathematical properties of the inverse 
k> ' and determinant of a large class of diagonally dominant matrices. 

5^ , By far, the most useful information about the inverses of such matrices in applications are bounds on 

their norms, so we discuss these properties first. A classical result of Varah [IB] states that if J is strictly 
diagonally dominant, i.e. if Aj(J) > for 1 < i < n, then the inverse of J satisfies the bound: 

|U -1 ||oo < max 



l<i<n Aj(J) ' 

Here, || • ||oo is the maximum absolute row sum of a matrix, which is the matrix norm induced by the infinity 
norm | ■ |oo on vectors in K n . Generalizations of this basic estimate can be found in [17], [H], and [9], but all 
involve the quantity maxi<i<„ 1/A,(J). In practice, however, one sometimes needs bounds when Aj(J) = 
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for some i, in which case the estimates appearing in |16i I17i 1121 [9] do not apply. A particularly interesting 
case is when A^(J) = for all i; we call such matrices diagonally balanced. 

We prove a tight bound on || J _1 ||oo for symmetric diagonally dominant J with positive entries that is 
independent of the quantities A* (J), and thus also of the maximum entry of J. Let S = (n — 2)I n + l n ln be 
the diagonally balanced matrix whose off-diagonal entries arc all equal to 1, and recall the Loewner partial 
ordering on symmetric matrices: A y B means that A — B is positive scmidcfinitc. We shall also write 
A>B it A — B is & nonnegative matrix. In Lemma 17.11 below, we show that if I > and J is a symmetric 
diagonally dominant matrix satisfying J > IS, then J y £S y 0; in particular, J is invertible. Throughout 
this paper, /„ and l n denote the n x n identity matrix and the n-dimensional column vector consisting of all 
ones, respectively. We also write I and 1 if the dimension n is understood. 

The following is our first main result. 

Theorem 1.1. Let n > 3. For any symmetric diagonally dominant matrix J with J^ > £ > 0, we have 

IIJ-'IU^IIS- 1 



£" "°° 2£(n - 2)(n - 1)' 

Moreover, equality is achieved if and only if J = £S . 

In other words, the map J i— > || J _1 ||oo over the (translated) cone of symmetric diagonally dominant 

matrices J > £S is maximized uniquely at the point of the cone; i.e., when J = £S. If the off-diagonal entries 

of J arc bounded above by m and the largest of the diagonal dominances Aj(J) is S, we also have the trivial 

lower bound: 

I < II /-ill 

2m(n-l) + S Sl1 l|o °' 

which follows from submultiplicativity of the matrix norm || ■ W^. Therefore, the value of || J _1 ||oc = 0(— ) is 
tightly constrained for bounded £, m and 5. 

To probe some of the difficulty of Theorem ! 1.1| we first derive an estimate from standard bounds in matrix 
analysis. The relation J y £S is equivalent to J -1 ■< (^S*) -1 = \S _1 [SJ Corollary 7.7.4], and therefore by a 
basic inequality [SJ p. 214, Ex. 14], we have || J _1 | < tHIS* -1 !! for any unitarily invariant matrix norm | ■ ||, 
such as the spectral, Frobenius, or Ky-Fan norms. It follows that 

l|.7-'IU<^l|.7-'lh<fl|S-'|b = ^Lj. (1) 

However, this bound is 0{— =), whereas the bound given in Theorem 11.11 is 0(— ). In some applications this 
difference can be crucial. For instance, we explain in Section [2] how Theorem 11.11 proves the consistency of 
the maximum likelihood estimator for some random graph distributions inspired by neuroscience |3J- 

As the following example shows, Theorem 11.11 fails to hold if we relax the assumption that J be symmetric. 

Example 1.2. For t > 0, consider the following diagonally balanced matrices and their inverses: 

(2 + t 1 1 + A 
Jt=\ 1 2 + t 1 + * 1 ; Jt 




r-l 1 



4 



t+3 


t-1 


1 + 1 


t+1 


t-1 


t+3 


t+1 


t+1 




Since some entries of J t 1 diverge, we have ||Jt 1 ||oo— > oo as < - >oo. □ 

One standard approach to proving norm estimates such as the one in Theorem 11.11 is a perturbation 
analysis. More specifically, given a symmetric diagonally dominant J with entries bounded below by £, one 
tries to replace each entry Jij by £ and prove that the norm of the inverse of the resulting matrix is larger. 
However, such a method will not succeed, even in the balanced case, as the following examples demonstrate. 



Example 1.3. Consider the following two balanced matrices: 



/12 4 1 7\ 
4 9 3 2 
13 7 3 


and H = 


/9 1 1 7\ 
16 3 2 
13 7 3 


^7 2 3 12J 




^7 2 3 12/ 



J = 



Here, H is obtained from J by changing its (1, 2)-entry to 1 and then keeping the resulting matrix balanced. 
A direct computation shows that ||P _1 ||oo < H^" 1 !!^- As another example, the following two matrices: 



J 




and H 




also have \\H 1 || oc < \\J 1 || o- Here, the (1, 2)-entry was changed without keeping the matrix balanced. □ 
As it will be an essential calculation for our proof of Theorem 11.11 we study (in Section |4j) the limit 

N= lim (S + tP)-\ (2) 

t— foo 

where P is a symmetric diagonally dominant matrix with Py <E {0, 1} and Aj(P) £ {0, 2}. Each such matrix 
P is a signless Laplacian of an undirected, unweighted graph G, possibly with self-loops. We compute the 
limit ([2]) in terms of the bipartition structure of the connected components of G. In particular, if G is 
connected and not bipartite, this limit N is the zero matrix (see Corollary 14. 5[) . For some recent work on the 
general eigenstructure of signless Laplacians, we refer the reader to [2] and the references therein. 
The following example reveals some of the combinatorics at work underlying Theorem 11.11 

Example 1.4. Consider the chain graph G with edges {1, n} and {i, i + 1} for i = 1, . . . , n — 1. If n is odd 

then N = since G is not bipartite, while if n is even, the limit N has alternating entries: 

v ~ n(n - 2) ' 
As another example, consider the star graph G, which has edges {1, i} for i — 2, . . . , n. In this case, 

1 ... 1 



N u = N zl 



2(n- l)(ra-2)' 



for i 



and Ni, 



2(n- l)(n-2) 



, otherwise. 



D 



We next describe our second collection of main results concerning the determinant and adjugate of positive, 
diagonally dominant symmetric matrices J. Recall the well-known formula: 



J- 1 = 



J* 



det( J) : 



in which J* is the adjugate (or classical adjoint) of J. To prove Theorem 11.11 one might hrst try to bound 
|| J*||oo from above and dct(J) from below, preferably separately. We first focus on estimating the latter, 
which is already nontrivial. 

It is classical that the determinant of a positive semidefinite matrix A is bounded above by the product 
of its diagonal entries: 



< dct(A) < Y[ Ai 



i=i 



This well-known result is sometimes called Hadamard's inequality [SJ Theorem 7.8.1]. A lower bound of this 
form, however, is not possible without additional assumptions. Surprisingly, there is such an inequality when 
J is diagonally dominant with positive entries; and, when J is also balanced, an improved upper bound. 



Theorem 1.5. Let n > 3, and let J be an nxn symmetric matrix with off-diagonal entries m > Jij > £ > 0. 
(a) If J is diagonally dominant, then 



n-\ 



det(J) / 1 fm /, m\\ ( 1 [m ( m 



nr=i ■& _ v 2( n - 2) v i \ i i) v 2 v t v * 

(b) // J is diagonally balanced, then 

det(J) / £ 2 \ 

nlu^- exp l^- 

The bounds above depend on the largest off-diagonal entry of J (in an essential way for (a); see Ex- 
ample [7T5]) , and thus are ill-adapted to prove Theorem ll.il For instance, combining Theorem 11.51 fa) with 
Hadamard's inequality applied to the positive definite J* in the obvious way gives estimates which are worse 
than (TT]). Nevertheless, Theorem 1 1 . 5 1 should be of independent interest, and we prove it in Section[7] 

As an immediate application of Theorems 11.11 and II. 5[ we obtain a bound on J* for balanced, positive J. 

Corollary 1.6. Let n > 3, and let J be a symmetric nxn matrix J with off-diagonal entries m > Jy > t > 0. 
If J is diagonally balanced, then the adjugate J* satisfies the bound: 

\\J*\\oo < 3n-4 c -l*/(4m 2 ) 



YltiJu ~ 2£(n-2)(n-l) 

Remark 1.7. We do not know how to prove this result other than by appealing to Theorems 11.11 and 11.51 

We finish this introduction with a brief overview of how our main results are proved. Theorem 11.11 will 
be generalized in Theorem 16.11 where we consider diagonally dominant matrices J > S(a,£) := al n +^1„1„ 
with a > (n — 2)£ > 0. We break up the proof of this general theorem into three main steps in Sections [3] 
to [6] where we write S(a, £) as S for simplicity. The first step considers the problem of maximizing || J _1 ||oo 
over symmetric diagonally dominant J with Jij and Aj(J) in some finite intervals (Section [3]). In this case, 
the maximum is achieved when J is on the corners of the space; namely, when J^ and Aj(J) are one of the 
endpoints of the finite intervals. In the second step (Section @|, wc analyze the behavior of corner matrices at 
infinity, showing that the limit of 11(5* + t-P) _1 ||oo as t — >• oo when P is a signless Laplacian (Pij £ {0, 1} and 
Aj(P) £ {0,2}) is at most ||5 _1 ||oo. Combined with the first step, this shows that the matrix S maximizes 
II t^ -1 II co over the space of symmetric diagonally dominant matrices J > S. The remainder of the argument 
deals with the behavior of || J _1 ||oo near S to show that S is indeed the unique maximizer (Section[5]). All 
three steps are combined in Section [51 The inequalities of Theorem 11.51 are proved using a block matrix 
factorization in Section [7J We conclude with a brief discussion of open questions in Section [5] 

2 Applications 

The (oo-norm) condition number Koo(A) = || -A||oo|| ^4~ 1 ||oo of a matrix A plays an important role in numerical 
linear algebra. For instance, the relative error \x — x|oo/|a;|oo of an approximate solution a; to a set of linear 
equations Ax = b is bounded by the product of k^A) and the relative size of the residual, \b — ^4i|oo/|^|oo 
(e.g., [5j p. 338]). Directly from Theorem ll.il we may bound the condition number of a positive, diagonally 
dominant symmetric matrix. Thus, numerical linear computation involving such matrices is well-behaved. 

Corollary 2.1. The condition number k 00 (A) of a positive, diagonally dominant symmetric nxn matrix A 
with largest off-diagonal entry m, smallest entry £, and largest diagonal dominance 5 satisfies: 

(2m(rc-l) + S)(3n-4) 
Kco{A > ~ 2£(n-2)(n-l) ' 



We next discuss another application of Theorem 1 1.1 1 to the numerical stability of inverses of a large family 
of functions. Let U be a convex open subset of R™, n > 3, and consider the map F: U — > R™ given by 

n 

F{x) = (F 1 {x),...,F n {x)), F i (x) = Y t g ij {xi + x j ), (3) 

3=1 

in which gij : R — > R are any continuously diffcrentiable functions with gij = gji, g\Ax% + Xj) > £ > 0, and 
g'ui^Xi) > for all 1 < i,j < n and x € U. 

Theorem 2.2. For any 9,9 e F(U), we have: 

■^ 3« — 4 ~ 

2«(n — l)(n — 2) 

Proof. The function i* 1 is continuously differentiable on [/, and its Jacobian J = Jf(x) satisfies: 

Jij(x) = ^ = 9ij{xi + Xj) > £, Jn{x) = * = 2g' u (2xi) + J^ g^fo + Xj) > ^ Jy (x). 



'.; 



o¥* 3& 



In particular, J is symmetric and diagonally dominant with off-diagonal entries bounded below by (. > 0. 
By the mean- value theorem for vector- valued functions p. 341], for each pair x,y € U, we can write: 

F(x) - F(y) = J(x - y), (4) 

in which the matrix J — L J(tx + (1 — t)y) dt is the element-wise average of the Jacobians on the line 
segment between x and y. Thus, J is also symmetric and diagonally dominant with Jy > i. In particular, 
J is invertible, which shows that the map F is invertible on U. By substituting 6 = F(x) and 6 = .F(y) into 
((4]), inverting J, and applying Theorem 1 1.1 1 to J -1 , we arrive at the inequality stated in the theorem. □ 

We close this section by explaining an application of Theorem 12.21 to probability and statistics. In [4], 
the maximum entropy distribution on (undirected) weighted graphs A = (Aij)fj =1 given an expected vertex 
degree sequence d = (di, . . . ,d n ) S R™ is studied, extending the work of Chatterjee, Diaconis, and Sly [1] 
in the case of unweighted graphs. When the graphs have edges in [0, oo), this distribution P# has random 
weights Aij that are independent exponential variables with E[Ay] = 1/(0* + 6j) > 0, where the parameters 
6 = (0i, . . . , 6 n ) £ R n are chosen so that the expected degree sequence of the distribution is equal to d. 

Given 9 with 1/m < ($i + 6j) 2 < l/£ for i j^ j, suppose we draw a sample graph G ~ Fg with weights 
A = (A^), and let di = J^j^i^ij be the degree sequence of G. The maximum likelihood estimator 6 for 6 
solves the moment-matching equations (the retina equations in theoretical neuroscience [3]): 

di = > — — — , for i = 1, ...,n. (5) 

i^ei+ej 

It can be shown that a solution 6 to (|5|) with di + 6j > 0, i =£ j, is unique almost surely and, more surprisingly, 
that the estimator 9 is consistent. In other words, the degree sequence of one sample from a maximum entropy 
distribution essentially determines the distribution for large n. We remark that solutions 9 to equations (0 
also exhibit a rich combinatorial structure which has been explored (within a more general framework) by 
Sanyal, Sturmfels, and Vinzant using matroid theory and algebraic geometry |llj . 

To see how Theorem 12.21 applies to this context, consider the function F determined as in ((3]) by setting 
gij(z) = — 1/z for i ^ j and ga = otherwise. In this case, equations ([5]) are concisely expressed as d — —F{&). 
Using Theorem 12.21 and arguments from large deviation theory, one can show that given any k > 1 and for 
sufficiently large n, we have the estimate: 



.. c-, 150\/m kloen , , , , ... 3 . . 

\8 - 0|co < j— \ ^— with probability > 1 t— r. (6) 



Thus, \9 — 6>|oo — > with high probability as n — > oo. Wc note that analogous arguments with the 0(-k=) 
bound from fl} fail to show consistency. We refer the reader to [4] for more details. 

3 Reduction to Corners 

To prove Theorem ll.il we need to show that the maximum of || J _1 ||oo over the space of symmetric diagonally 
dominant matrices J > S := al n + £l„lj l is achieved at J = S. A priori, it is not even clear that a maximum 
exists since this space is not compact. In this section we consider maximizing | J _1 ||oo over compact sets of 
symmetric diagonally dominant matrices J, and we show that the maxima occur at the corners of the space. 
In subsequent sections we analyze these corner matrices in more detail. 

Fix m > 1, and let T> = T> m denote the set of n x n matrices of the form J = S + (to — £)P where P is 
some symmetric diagonally dominant matrix satisfying 

< Pij < 1 for i^j and < A 4 (P) < 2 for i = 1, . . . ,n. 

We say that J G T> is a corner matrix if 

Pij E {0,1} for i^j and A.;(P)e{0,2} for i = l,...,n. 

Equivalently, J is a corner matrix if P is a signless Laplacian matrix. Let T denote the set of matrices J £ T> 
that maximize ||J ||oo- This set is closed and nonempty since the function J n- || J _1 ||oo is continuous and 
T> is compact in the usual topology. Let ei, . . . , e n be the standard column basis for R n , and set en = eiej 
and e.ij = (e-i + ej)(e, + ej) T for all i ^ j. Our main result in this section is the following. 

Proposition 3.1. Every J € T is path- connected to a corner matrix. 

Proof. Let J e T. Wc will show that if £ < Jij < m for some i ^ j, there is a path in T from J to a matrix 
J' that differs from J only in the (i, j)-entry, with J[- G {£, m}. Similarly, if < Aj(J) < 2(m — £) for some 
1 < i < n, then we find a suitable J' with A^J') € {0, 2(to — £)} that differs from J in the (i,i)-entry. 
Repeatedly applying these steps, it follows that there is a path in T from any J € T to a corner matrix. 
For the first part, suppose that £ < Jy < m for some i ^ j. Consider the nonempty, closed set 

W = {J + teij :£<J ij +t<m} n T. 

We claim that W contains a matrix J' with J|- G {^, m}. Suppose not, and let J' e W be a matrix with 
minimum (i, j)-entry. By Proposition 13.21 (a) below, J' + fey 6 T for all t in a small neighborhood of the 
origin. Thus, there is another matrix in T that has a smaller (i, j)-entry than J', a contradiction. The proof 
for the other part is similar (using Proposition 13.21 (h)). □ 



To complete the proof of Proposition 13. 11 it remains to show the following. 
Lemma 3.2. Let J 6 T and i ^= j be distinct indices in {1, . . . , n}. 

(a) If 1 < Jij < m, then J + te^ G T for all t G M. in some neighborhood of the origin. 

(b) If < Ai(J) < 2(to — 1), then J + tea € T for all t G R in some neighborhood of the origin. 

Proof. We only prove (a) as (b) is analogous. Suppose that 1 < Jy < m for some t ^ j. Let A' = J^ 1 and 

set Ai, . . . , K n to be its columns. Also, let q e argmaxi< p < n | A p |i and I = {1 < p < n: A p( j 7^ 0} so that 

5^ \K Pq \ = |A g | 1 = ||A|| 00 = II J- 1 ]^. 

pel 

By the Sherman-Morrison- Woodbury formula, we have 



where rjij := (e,; + ej) T K{e.i + ej) > since K y 0. The formula implies that for sufficiently small e > and 
all |*| < e, the (p, g)-entry of ( J + iey) -1 has the same sign as K pq ^ for all p £ 1. Let us further suppose 
that e is small enough so that J + ie.y e 2? and 1 + rjijt > for all |*| < e. The 1-norm of the q-th column 
(J + teij)^ 1 of the matrix of (J + ie.^-)" 1 now satisfies: 



(</ + teij)- 1 ^ = || J- 1 !^ - ----- + --A, 

J7? ' " " 1 + Vijt 1 + Vijt 



for all \t\ < e 



(K iq + Kj q ) Y, p ei s ^ n ( K pq)(. K ip + K w) and ^ := l A ^<z + ^'«l S P ^i l-^ip + ^j'pI- Recalling 



where 0$. 

that J £ T achieves the maximum || J _1 ||oo, it follows that 



This inequality implies that ipij = 4>ij = 0, so ||(J + te 



'i j i 



= 11 J- 



for all |i| < e as required. 



□ 



4 Corner Combinatorics and Exact Limit Formulae 

In the previous section, we saw that corner matrices played an important role in the optimization of || J _1 ||oo 
over all matrices J > S := al n + ^1„1„ . A corner matrix may be expressed as: 

J = S+ (m - i)P, 

in which P is a symmetric diagonally dominant matrix with Pjj 6 {0, 1} for all i ^ j and Aj(P) £ {0,2} 
for 1 < i < n. Every such matrix P is a signless Laplacian of an undirected unweighted graph (possibly 
with self-loops) G = (V,E) on vertices V = {1, . . . ,n} and edges £?. That is, we can write P = D + A 
where D is the diagonal degree matrix of G and A is its adjacency matrix. We study the limit ([2]) using the 
combinatorics of the graphs associated to the matrix P. 

Example 4.1. Let S — (n — 2)I n + l n l^ and G be the chain graph from Example 11.41 For n = 4, 



S + tP 



/3 + 2* 1 + * 

1 + * 3 + 2* 

1 1 + t 
\1 + * 1 



1 l + t\ 

1 + t 1 

3 + 2* 1 + * 

1 + * 3 + 2*/ 



with inverse: 



(S + tP)- 1 = 



1 



4(2* + 3)(* + l) 



/* 2 + 5* + 5 -(i+1) 2 i 2 +i-l 

-(* + l) 2 t 2 + 5* + 5 -(* + l) 2 

* 2 +*-l -(* + l) 2 * 2 + 5* + 5 

\-(*+l) 2 * 2 +*-l -(* + l) 2 



-(* + l) 2 \ 
t 2 + * - 1 

-(i + 1) 2 
t 2 + 5* + 5/ 



Each entry of (5 + tP)~ l is a rational function of t, with numerator and denominator both quadratic. Thus, 
each entry converges to a constant as * — > oo, and from the expression above, we see that the limit matrix N 
has entries Nij = (— l) t+J /8, as predicted by the formula in Example 11.41 
If one adds the edge {1, 3}, the corresponding matrix (S + tP)^ 1 is 



(s + tpy 



4(i + 3)(* + l) 
Thus N — 0, as the graph is no longer bipartite (see Corollary 14. 51 below) 



/2t + 5 -t-1 -1 -t-l\ 

-* - 1 3* + 5 -t-1 t-1 

-1 -t-1 2* + 5 -t-1 

\-t-l t-1 -t-1 3* + 5/ 



□ 



We begin with the following simple fact that allows us to invert certain classes of matrices explicitly. 

Lemma 4.2. The following identity holds for any a =/= and £ ^= —a/n: 

1 £ 

(al n + ll n l ) = l n - ■ - -1„1 . 

a a{a + In) 

Given a signless Laplacian P, let L € ]R nx l £ l be the incidence matrix of the graph G associated to P; 
that is, for every vertex v £ V and edge e € E, we have: 

1 if v is in e and e is not a self-loop, 

V2 if v is in e and e is the self- loop (v, v), 
^ otherwise. 

Consequently, P = LL T . Using this decomposition of P, we derive the following formula for N. 
Proposition 4.3. The limit N in (0) satisfies: 

N = S- l -S- l ' 2 {X T yx T S- 1 l\ 

where X = S~ 1 > 2 L and (X T Y is the Moore-Penrose pseudoinverse of X T . Furthermore, NL = 0. 

Proof. Using the Sherman-Morrison- Woodbury matrix formula to expand (S + tLL )~ 1 , we calculate: 

iV = limt-xx, (S + tLL T )~ 1 

= lim^^ [ S- 1 - S- l L(t- l I + L T S- 1 L)- 1 L T S- 1 } 
= S- 1 - S- 1 ' 2 [lim^eo X(t- l I + X T X)- 1 } X T S- 1 ' 2 

= s- 1 -s- 1 ' 2 (x T yx T s~ 1 ' 2 , 

where in the last step we used an elementary pseudoinverse identity in matrix analysis [3 p. 422, Ex. 9]. To 
show that NL = 0, we note that N is symmetric and that 

[NL) T = L T N = X T S~ 1/2 - X T (X T )^X T S~ 1/2 =0. □ 

Let Ni, . . . , N n denote the columns of the matrix TV. The following is immediate from NL = 0. 

Corollary 4.4. For each edge {«, j} of G, we have Ni = —Nj. In particular, Ni = for each self-loop {i, i}. 

Suppose that the graph G has connected components G\, . . . , G&. After relabeling the vertices, both P and 
L are block-diagonal with matrices Pi, . . . , P& and L\, . . . , Lk along the diagonal. Furthermore, Pj = LiLj 
for each 1 < i < k. The components of G also induce a block-structure on the limit N, and we denote these 
blocks by N[i,j] for i,j £ {1, . . . , k}. These blocks display many interesting symmetries. Firstly, the entries 
in each block are all equal up to sign. Secondly, the signs in a block N[i,j] depend on the bipartite structures 
of the components Gi and Gj. We say that a bipartite graph is (jp, q)-bipartite if the partitions are of sizes p 
and q, respectively. Note that bipartite graphs cannot have self-loops. 

Corollary 4.5. In each block N[i,j], the entries are equal up to sign. If the i-th component Gi of G is not 
bipartite, then N[i, j] = N[j, i] = for all j = 1, . . . , k. In particular, N = if and only if all components 
Gi are not bipartite. Suppose Gi is (pi,qi) -bipartite and Gj is (j?j,qj) -bipartite. Then, after relabeling the 
vertices, the matrix N[i,j] has the block structure 

/ 1 1 T -1 1 T \ 

I Pi p Pi Q \ 

N[i,j] = Cij \ ' T 3 , for some constant Cij £ R. 

\ ^ 1 1. 1 p j 1 9. 1 ? J / 



Now that we understand the block structure of N, we want to compute the constants c.y ■. Our approach 
is to simplify the formula in Proposition 14.31 by expressing the incidence matrix L in a more suitable form. 



Proposition 4.6. Let G be a connected graph with n vertices. If G is not bipartite, then rankL = n. 

Proof. It suffices to prove that the positive semidefinite matrix P = LL T has rank n; i.e., if x £ K™ satisfies 
x T Px = 0, then x = 0. Write 



X PX — 2_^ 2_^, "ij X i X j ~ 2-~t 

{i,j}eE 



so that x T Px = implies 



^ i 



if {i, i}6£ and 



(x, l +x j f+ Y, 2x l 

i: {i,i}£E 



i£{i,j] £E. 



(7) 



If G has a self-loop {i, i}, then Xj = by 0. If G has no self-loops, there is an odd cycle (ii, «2, ■ • ■ , *2m+i) 
for some m > 1 (because G is not bipartite). Applying condition to each successive edge in this cycle 



shows Xi x 



V l2m + 1 



0. In either case Xi — for some vertex i £ V. A repeated application of Q 



then reveals that Xj = for all vertices j connected to i. Since G is connected, we have x — as desired. □ 
Recall that e, € K." denotes the «-th standard basis column vector. 

Proposition 4.7. Lei G be a connected bipartite graph on n vertices. Then rankL = n — 1. Lei [/ 6e i/ie 
n x (n — 1) matrix with columns e.\ + cr-iei /or 2 < i < n, where o~i = — 1 i/ vertex i is in the same partition 
as vertex 1 and u, = 1 otherwise. Then, L = C/L? /or some (n — 1) X |L7| matrix B of rank n — 1. 

Proof. Recall that the columns of L are of the form e, + e^ where {«, j} € L7 (there are no self-loops because 
G is bipartite). There is a path (1, i±, i<i, ■ • • , i m ,j) from vertex 1 to vertex j for each 2 < j < n, so 

ei + (-l) m ej = (ei + e^J - (e n + e l2 ) + • ■ ■ + (-I) m (e lm + e^). 

Thus, rankL > n — 1. Conversely, for each edge {i, j} where vertex i is in the same partition as vertex 1, 

e, + ej = -(ex - e,) + (e 1 + ej), 



so rank L < n — 1. This equation also allows us to write L 
on B follows from that of L. 



UB for some matrix B, and the rank condition 

□ 



Recall that G has k components Gi, . . . , Gk and L is block-diagonal with matrices Li, . . . , Lfc. If Gi is 
bipartite, we write Li = UiBi as in Proposition ^. 71 If Gi is not bipartite, we write Li = UiBi where Ui = I is 
the identity matrix and B, t = Li. Let r be the number of components of G which arc bipartite. If U and B arc 
block-diagonal matrices constructed from U\, . . . ,Uk and B\, . . . , B^, then L = UB where U £ ]R rax (™ _r ) and 
B £ M( ,l ~ r ) x l £ l both have rank (n — r). Note that U contains information about the sizes of the bipartitions 
of each component whereas B contains information about the edges. Let U< p „\ denote the matrix 



U, 



(p,q) '■' 



After relabeling the vertices, we have 



U 



( U (pi,qi) 




U, 



(P2.92) 






—Ip-1 





1 





u. 



(p r ,q r ) 



\ 






I. ) 



u 







I, 



(8) 



where s = n — $^i=i(P» + ?i) is the total number of vertices in the non-bipartite components of G. Our next 
result shows that dependence on the matrix B can be removed in Proposition 14.31 This new formula for N 
also gives us a method to compute its entries explicitly. 

Proposition 4.8. Let U be as in {Sp. The limit N in Hj] satisfies: 

n = s~ 1 - s- 1 u(u T s~ 1 uy 1 u T S' 1 . 

Therefore, N depends only on the sizes of the bipartitions of each component Gi . 

Proof. Write L = UB and P = U(BB T U T ). First, note that BB T is positive definite since B € R("- r ) x l B l 
has rank (n — r). Now, U T S~ 1 U is positive semidefinite, being a congruence of a positive definite matrix S^ 1 
(Sylvester's Law of Inertia). Moreover, since U is injective as a linear map and S 1-1 y 0, for any x <G M. n ~ r : 



x T U T S- 1 Ux = 



Ux = 



0. 



Thus, U T S~ 1 U is positive definite. Then the eigenvalues of BB T U T S^ 1 !/ are all positive, being a product 
of two positive definite matrices (see [71 Lemma 2]). Therefore, the matrix t~ l I + BB T U T S~ 1 U is invcrtiblc 
for all t > 0, and so by the Sherman-Morrison- Woodbury formula, we have 

(S + tP)- 1 = (S + tU{BB T U T ))- 1 

= S- 1 - S^Uff- 1 ! + BB T U T S- 1 U)- 1 BB T U T S-\ 



Taking limits of this equation as t — > oo, the result follows. 
To state our ultimate formula for N, let us first define: 



□ 
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E 



Pi + q% 



Pi +<li 



(l^.-lj) 



— (1 T -1 T ) 0-l T 



Y = 



1 



' 1 1 T 
*-pi ± pi 

Pi +11 \ -IgilJ 



V 



-1 1 T 
1 1 T 

± qi ± q 1 



1 1 T 
-1 1 T 












Proposition 4.9. Set"/, y, andY as above, which depend only on the bipartite structures of the components 
of the underlying graph G. We have the following formula for the limit in (0): 



N 



lim (S + tP)- 

t—toc 



-Y 



t 



;yy 



a(a + £7)' 

Proof. We outline the computation of N. For simplicity, let us write S 1-1 = al n 

U(U T S- 1 U)- 1 U T = U(U T (aI n - bl„ll)U)- 1 U T 
= U(aW -bvv T )- l U T , 



bl n l±. Then, 



(9) 



10 



where W = U T U and v = U T l n - By the Sherman- Morrison- Woodbury identity, we have that 

(aW-bvv T )- 1 =a- 1 W- 



„-i w -i a^W-H-bv^a^W- 1 



l + v T a- 1 W- 1 (-bv) 

= -W- 1 + — -—W- 1 vv T W^\ 

a a(a — bq) 

where <r = v T W~ 1 v (it is easy to show that a — bq > 0). Substituting this back into (|9|) gives us 

U(U T S- 1 U)- 1 U T = -Z+ . b , , zz T , (10) 

a a{a — o?) 

where Z = U(U T U)~ 1 U T , z = Zl n , and q = l^Zl n . From the block-diagonal structure of U, the matrix Z 
is also block-diagonal with blocks 

U^UJU^UJ = I Pi+qi - -i_ ( lp ^T " 1P ;?V *»i = l,...,r, 

Pi-t qi y—i qi i- p . *-qi*- qi J 

where we have computed {UjU^ 1 using Lemma T4. 2 1 Next, one shows that Z = I n — Y, z — l n — y, and 
q = n — 7. Finally, substituting 

1 ^ 

S^ 1 = al„ -bl n lZ, a = -, 6=— ? — — — -, (11) 

a a(a + in) 

and equation (JTOJl into Proposition 14.81 gives the desired result. D 

Corollary 4.10. For 1 < i,j < r, the constants Cij in the blocks N[i,j] in Corollary \4-.5\ are: 



a , (pi - qi) 



T+7 



a(a + h){Pi + <li) \ £ Pi + Vi 

f Pi-qi\ (pj-Qj 



-^ a {a + i 1 )\p i + q i )\p j +q j i " 

Finally, we write down an explicit formula for H-ZVlJoo and verify that ||iV||oo < US' - ||oo- 
Corollary 4.11. Ifr = 0, then ||JV||oo =0. Ifr>l, let d = £T =1 \Pi - Qi\- Then 

a aya + £7) i<t<r p^ + q t 

Proof. Recall that if r = (i.e. no component of G is bipartite), then N = by Corollary |4.5l Now if r > 1, 
we may assume that pi > qi for all i after a relabeling of vertices. Observe that ca > and Cy < for all i 
and j 7^ i. Indeed, this follows from a > 0, £ > 0, and 

7 = y- fa-%) 2 >p (p«-g') 2 = y- (Pj-gj) 2 > 

Consequently, the 1-norm of rows in the i-th block of N is 

f , x v-, , x 1 l(pi-g i )(rf-2p i + 2g i ) 

(ft + qi)Cii - 2j<Pj + Qj)Cij = - + -. — —fx, ; — \ — . ( 12 J 

^ a a{a + £-f){p l + qi) 

and H-ZVlloo is the maximum of these 1-norms. D 
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Corollary 4.12. For all signless Laplacians P, we have 



N 



\im(S + tPy 



<\\s- 



1 , lfo-2) 
a a(a + In) ' 



with equality if and only if P is the zero matrix. 

Proof. Let N = lim f __ >00 (S' + tP)^ 1 . If r = then ||iV||oo = 0, so suppose that r > 1. As before, assume that 
Pi > Qi for all £ = 1, . . . , r. It suffices to show that the 1-norms of the rows of N, as computed in (fT2|) . are at 
most ||5 _1 |joo7 with equality achieved only at J = S. The inequality is trivial ii Pi = qi so we may assume 
Pi — Ii > 1- We outline the proof and leave the details to the reader. The key is to show that 



{a + £n)(d-2pi + 2qi) 



P, 



< (a + £n)(d-2) < (a + £j)(n-2). 



The latter inequality is equivalent to 

< {n-d){a + 2l-ln) + 2l{d-7) + nl(n-2d + ~f). 
The first summand is nonnegative because S is diagonally dominant, while the last summand satisfies 



n - 2d + 7 



y-^- > o. 



Finally, if equality H^VJIoo = 115 1 \\ OQ is achieved, then s = and qj = for all j, so P = 0. 



D 



5 Analysis of || J l^ in a Neighborhood of S 

The arguments in Sections [3] and U show that for J > S :— al n +il„l„ , the maximum of || J _1 j|oc is attained 
at J — S. To prove that S is the unique maximizer, we will show that ||«/ _1 ||oo is strictly decreasing near S. 
Let P > be a nonzero symmetric diagonally dominant matrix, and consider the function 



f(t) = \\(S + tP)^\\ 



t> 0. 



In our proof, we study the linear part S l — tS 1 PS x of the Neumann series for (S + tP) 1 . Let us define 
g(t) = WS^ 1 — tS~ 1 PS~ 1 1 1 oo and h(t) = f(t) — g(t). Our main result in this section is the following. 

Proposition 5.1. The function f{t) is differ entiable at t — and /'(0) < 0. 

Proof. Since f(t) — g(t) + h(t), the result follows from Propositions 15.21 and 15.31 □ 

Proposition 5.2. The function h{t) is differ entiable at t — and h'(0) = . 

Proof. For sufficiently small t > 0, by the Neumann series for (I + tS~ 1/>2 PS~ 1 / 2 )~ 1 we can write 

(s + tpy 1 = s- 1 / 2 ^ + ts- 1 ' 2 ps- 1 / 2 r 1 s- 1 i 2 = s- 1 ' 2 (JTi-tfis^^ps- 1 / 2 )"] s^ 1 ' 2 . 

\k=0 / 

By the reverse triangle inequality and submultiplicativity of || • ||oo, we have 

\h(t)\ < ||(5 + tP)- 1 - S- 1 + tS^PS^W^ 

—t)"{S -'-I'd -'- 1 

\k=2 
£-1/2 (5-l/2p 5 -l/2) 2 gl/2 (5 + tp) -l 



s- 1 ' 2 yi-tfis-'^ps- 1 / 2 )" S" 1 



/2 



<2i 2 ||5- 1 P5- 1 P|| co ||5- 1 | 
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where the last inequality holds for sufficiently small t > since by continuity \\(S + tP) noo ^ *\\u noo 
small t. It follows that h'(0) = lim t ^o(h(t) - h(0))/t = 0. □ 

Proposition 5.3. The function g(t) is differentiable at t — and g'(0) < 0. 

Proof. Set Q = S~ 1 PS~ 1 . Note that for sufficiently small t > 0, the entries of 5 1-1 — tQ have the same sign 
as the corresponding entries of S _1 . Since (S~ l )u > and (5 ,_1 )y < for i ^ j, we can write 

g(i) = ||5 _1 - tQHoo = max,; Ej \( s ~ 1 )i] ~ f Qij\ 

= max, [ (S- 1 ^ - tQ u + E, # , (tQ w - (5- 1 ) ii ) ] 

= max, [ - (Q« - X) i#i Q« ) * + ll^ _1 ||oo] 

= -^+||5- 1 ||oo. 

where £ = mini (Q« — Ej^j Q»j) > by Proposition 15.41 below. Thus, g'(0) = — £ < as required. □ 

Proposition 5.4. Lei Q = S^PS' 1 . Then Q u - E, ¥ ; Qy > /or all i. 
Proof. For simplicity, let us write S~ l = al n — 61„1„. Then, 

Q = {al - 611 T ) P (a/ - 611 T ) = a 2 P - abpl T - ablp T + 6 2 7rll T 
where p = PI = (pi, . . . ,p n ) and tt = 1 T P1. It is straightforward to check that 

Qu - Ej^i <2y = a 2 ( 2 -P« -Pi) + ab P i{n - 4) + 67r(a + 26 - 6n). 
From equation (|11[) . we get a/6 = a/£ + n. Substituting this relation and rearranging gives us 

Qu - E^ Q« = ^ [( j + n 2) (| + 4) + 4] A,(P) + 






2 (f + .-!)(„- 3) + f]fl, + *(f + j) Ef,, 



j,k=£i 



Because S 1 is diagonally dominant, we have a/i > n — 2 > 0. It is not difficult to deduce that if P 7^ 0, then 
the above expression is always positive, as required. □ 

6 Proof of Theorem 11.11 

Theorem 1 1.1 1 is a special case of the following theorem when a = (n — 2)1. 

Theorem 6.1. Let n > 3 and suppose S = al n + ^l n l^ is diagonally dominant with a, £ > 0. For all n x n 
symmetric diagonally dominant matrices J > S , we have 



|J _1 |loo < Halloo = 



a + 2£(n-l) 



a (a + in) 

Furthermore, equality is achieved if and only if J = S . 

Proof. Recall from Section[3]that T> m is the set of symmetric diagonally dominant matrices J with I < Jij < m 
and Aj(S) < Aj( J) < Aj(5) + 2(m — (.) for all i ^ j. Recall also that S + (m — t)P G V m is a corner matrix 
if P is a signless Laplacian. Let T m be the set of matrices J e T> m maximizing ||J _1 ||oo- We claim that for 
sufficiently large m > £, we have 7m = {S}. Indeed, from Corollary 14. 121 and for large m: 

||^ -1 ||oo < !|>5 1 ||oo, for ah corner matrices J <G V m \ {S}. 

Thus, by Proposition l3.11 every J £ T m must be path-connected to the corner matrix S. Since Proposition l5.ll 
implies that S is an isolated point in T m , we must have T m = {S} as claimed. 

Finally, suppose J* > S is a symmetric diagonally dominant matrix with ||(J*) _1 ||oo > ||<S _1 ||oo- We will 
show that J* = S, which proves Theorem 11.11 We assume that m is sufficiently large with J* € T> m and 
Tm = {S}. Then S is the unique maximizer of || J _1 ||oo for J G 2?m so that J* = 5, as desired. D 
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7 Extensions of Hadamard's Inequality 

Our arguments for proving Theorem [T3] are inspired by block LU factorization ideas in J3J. For 1 < i < n, let 
Ju\ be the lower right (n — i + l)x(n — i+1) block of J, so Jm = J and J(„) = (J rm ). Also, for 1 < i < n — 1, 
let 6(j) G M™~ 4 be the column vector such that 

J (i) = ( Jil b ^ 

Then our block decomposition takes the form, for 1 < i < n — 1 , 

1 U {1) \ ( Si 

with 



J w lo j„l'J Um J, 



« ^(i+i) 



r A 6 « J (i+i) 6 w\ , 

Jii 1 V - and U (i \ = b h\J,. 



; « -°(i) J (%+!)• 



Notice that det(J) = JnnY[i=i s ^ or equivalently, 



n— 1 re— 1 / f,T t — 1 l 

n£=n i-^|. as) 



dg(J) = TT il = TT ( ! _ TO)°M 

riz=i j « i=1 j " i=1 v 



It remains to bound each factor Su/Ju. We first establish the following results. 

Lemma 7.1. Let J be a symmetric diagonally balanced n x n matrix with < £ < Jy < m /or i =£ j. Then 
£S X J X mS, and i/ie eigenvalues Xi < ■ ■ ■ < X n of J satisfy 

(n - 2)£ < A s < (n - 2)m /or 1 < i < n - 1 anrf 2(n - 1)^ < A„ < 2(n - l)m. 

Moreover, if J is diagonally dominant, then the lower bounds still hold. 

Proof. We first show that if P > is a symmetric diagonally dominant matrix, then P >: 0. For any x € R", 



.r _ 



p.x = £ p^ 2 + 2 ^ p ijXiXj >J2[J2 p ij) x i +2 Yl p v x ^ = E fate + x of ^ °- 



Since the matrices P = J — £S and Q = mS — J are symmetric and diagonally balanced with nonnegative 
entries, it follows that P, Q >z by the discussion above, which means £S < J < mS. The eigenvalues of S 
are {n — 2, . . . ,n — 2,2(n — l)},so the result follows by an application of [5] Corollary 7.7.4]. If J is diagonally 
dominant, then £S ;< J, and hence the lower bounds, still holds. □ 

Lemma 7.2. Let J be a symmetric diagonally balanced n x n matrix with < £ < Jij < m for i ^ j . For 
each 1 < i < n, let J^ be the lower right (n — i + l)x(n — i + 1) block of J as defined above, and suppose the 
eigenvalues of Ju\ are Ai < • • • < \ n -i+\ ■ Then 

(n — 2)£ < Xj < (n — 2)m for 1 < j < n — i and (2n — i — 1)£ < A n _.; + i < (2n — i — l)m. 

Moreover, if J is diagonally dominant, then the lower bounds still hold. 

Proof. Write J(j) = H + D, where H is the (n — i + l)x(n — i + l) diagonally balanced matrix and D is diagonal 
with nonnegative entries. Note that (i — 1)£I ;< D ^ (i — l)ml, so (i — 1)£I + H < J^ ^ (i — l)ml + H . 
Thus by [5] Corollary 7.7.4] and by applying Lemma 17. II to H, we get, for 1 < j < n — i, 

(n -2)£={n-i- l)£ + (i - l)£ <Xj<{n-i- l)m + (i- l)m = (n - 2)m, 
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and for j = n — i + 1 , 

(2n - i - 1)£ = 2(n - i)£ + (i - 1)£ < A„_ l+ i < 2(n - i)m + (i- l)m = (2?i - i - l)m. 

If J is diagonally dominant, then (i — \)£I + H < Ju\ and hence the lower bounds still hold. □ 

Proof of Theorem 1 1.5t For part (a), suppose J is diagonally dominant. For each 1 < i < n — lwe have 
J a > yi,-4j Jij > ^Z\ln-ii an d by Lemma 17.21 the maximum eigenvalue of J,~ ., is at most , \\e - Thus, 



^J,M>iw 1 blsbu\ i &Zs 6r,^ ./fn-i + llm \/ 6 r,/i & 



3 «<Vi) 6 « _^ 6 (*) & W 1 & w 6 (») v/(n - j + 1) m V (0 (l) 



Jii "(n-2)^ J, ~ (n-2)£ bj^l ~ (n - 2)1 bj^l 

Since each entry of bu\ is bounded by £ and to, the reverse Cauchy-Schwarz inequality [151 Ch. 5] gives us 



\i) J (i+i)"(i) y/{n -i + 1) m £ + m _ 1 [rn f m 



Jii {n - 2)£ 2y/em(n-i + l) 2(n - 2) \ £ \ £ 

Substituting this inequality into (|13p gives us the desired bound for part (a). 

For part (b), suppose J is diagonally balanced, so Ju < (n — l)m for each 1 < i < n — 1. By Lemma 17^1 
we know that the minimum eigenvalue of J,~ ... is at least (2n-l-2)m > so 

6 (i) J (7+i) 6 W > 1 & W 6 (') > 1 (n-z + l)£ 2 _ (n-i + 1)^ 2 



Jii (2n — i — 2)to Ji 2(n — 1)to (n — l)m 2(n — 1) 2 to 2 

Substituting this into (|13p and using the inequality 1 + a; < cxp(x) gives us the desired bound: 

dct(J) ^r 1 / {n-i + l)£ 2 \ ( {n + 2)£ 2 \ ( £ 2 \ 

Yf^j- * g ex P h^DwJ = CXP (^(n-ljm'j * CXP [-4^)- D 

We close this section with several examples. 
Example 7.3. The matrix S — (n — 2)I n + l n l^ has eigenvalues {n — 2, . . . , n — 2, 2(n — 1)}, so 

det(S) 2(n-2) n - 1 (n-l) ( 1 \" _1 2 



Example 7.4. When J is strictly diagonally dominant, the ratio det(J)/J|™ =1 Ju can be arbitrarily close 
to 1. For instance, consider J = al n + lnl„ with a > n — 2, which has eigenvalues {(n + a), a, . . . , a} so 

dct(J) (n + a)^™" 1 _ 

1 1 as a -> oo. D 



nr=i ■& " (« + 1)" 

Example 7.5. The following example demonstrates that we require an upper bound on the entries of J in 
Theorem [T3] (a). Let n — 2k for some k £N, and consider the matrix J in the following block form: 

J =(i T)> A = (km + k£ - 2£)I k + £l k lj, B = ml k lJ. 



B A / 

By the determinant block formula (since A and B commute), we have 
det(J) = det(A 2 - B 2 ) 



det 



(km + k£- 2£) 2 I k + (2k£m + 3k£ 2 - 4£ 2 - km 2 )l k ll 



4£(k - l)(km + k£-£)- (km + k£ - 2£) 2k - 2 1 
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where the last equality is obtained by considering the eigenvalues of A 2 — B 2 . Then 

det(J) 4£(k - l)(km + k£ - £) ■ (km + k£ - 2£) 2k ~ 2 4£ ( 2£ \ 

1 _ rra = — — — - — ttt > cxp — as k ->• oo. 

Ui=iJu (km + ki-i) £ + m V t + mj 

Note that the last quantity above tends to as m/t — >• oo. D 

8 Open Problems 

As an analogue to Theorem 16. 1[ we also conjecture a tight lower bound for symmetric diagonally dominant 
matrices J > whose off-diagonal entries and diagonal dominances are bounded above. Observe that when 
Jij < tn and Aj(J) < 6 for all i ^ j, then J < (m(n — 2) + S)I n + ml„l„ . 

Conjecture 8.1. Let n > 3 and let S(a, m) = al n + ml n l„ . For all n x n symmetric diagonally dominant 
matrices < J < S(a,m), we have 

n r— in ^ncv n— in a + 2m(n-l) 

|| J ||oo > \\S(a,m) || 



oo 



a (a + mn) 

Moreover, equality is achieved if and only if J = S(a,m). 

In working towards the proof of this conjecture, the following problem may be useful. 

Problem 8.2. Given a signless Laplacian P of a graph G, give an exact combinatorial formula for the entries 
of (S + tP)~ 1 for any t > 0. More precisely, since each entry of (S + tP)~ 1 is a rational function oft, derive 
a formula for the coefficients of this rational function in terms of the combinatorics of the graph G. 

We also conjecture that a dependence on the largest entry can be removed in Theorem II .51 (b). 

Conjecture 8.3. For a positive, diagonally balanced symmetric J , we have the bound: 

det(J) < dct(S) _ 2 A 1 V" 1 k 2 



rr=i J « " (n-l)" V n-lj e 
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